On the supernova remnants with flat radio spectra 
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Abstract A considerable fraction of Galactic super- 
nova remnants (SNRs) characterize flat spectral indices 
(a < 0.5). There are several explanations of the flat 
radio spectra of SNRs in the present literature. The 
most of models involve a significant contribution of the 
second-order Fermi mechanism but some of them also 
discuss high compressions (> 4), contribution of sec- 
ondary electrons left over from the decay of charged 
pions, as well as the possibility of thermal contamina- 
tion. In the case of expansion in high density environ- 
ment, intrinsic thermal bremsstrahlung could theoret- 
ically shape the radio spectrum of an SNR and also 
account for observable curved - "concave up" radio 
spectra of some Galactic SNRs. This model could also 
shed a light on the question of fiat spectral indices de- 
termined in some Galactic SNRs. On the other hand, 
present knowledge of the radio continuum spectra (in- 
tegrated flux densities at different frequencies) of SNRs 
prevent definite conclusions about the significance of 
proposed models so the question on flat spectral in- 
dices still remains open. New observations, especially 
at high radio continuum frequencies, are expected to 
solve these questions in the near future. Finally, as 
there is a significant connection between the majority 
of Galactic SNRs with flat integrated radio spectrum 
and their detection in 7-rays as well as detection of ra- 
diative recombination continua in their X-ray spectra, 
the analysis of high energy properties of these SNRs is 
very important. 
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1 Introduction 

The spectra of SNRs in radio continuum are usually 
represented by a power-law, reflecting the pure syn- 
chrotron radiation from the SNR shell: 



S^cnv 



(1) 



where Si, is the spatially-integrated flux density and a 
is the radio spectral index. 

Test-particle diffusive shock acceleration (DSA) the- 
ory predicts that for strong shocks, radio spectral index 
a is approximately 0.5 (Bell 1978ab). Indeed, mean 
value of the observed radio spectral indices, for Galac- 
tic SNRs, is around 0.5 and corre spondent distribution 
is roughly Gaussian (jGreenl 120091 ). On the other hand, 
the dispersion in that distribution is non-negligible (see 
Fig. 1 in Reynolds et al. 2012). 

In the case of shocks with sufficiently low Mach num- 
ber (less than around ten), steeper spectra (a > 0.5) are 
generally expected. In test-particle DSA, for the syn- 
chrotron spectral index a, the following relations hold: 
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where 7 is the post-shock thermal gas adiabatic index, 
X is the compression ratio for parallel shock, M is the 
(upstream adiabatic) Mach number that represents the 
ratio of shock velocity to the local sound speed and s 
is the energy spectral index. On the other hand, only 
few old SNRs would b e expected to have such weak 
shocks (JRevnoldsluOllI) . In fact, SNRs wit h steep spec- 
trum (g > 0.5) are usually youn g objects (|Greenll2009 : 
Gao et al.ll201lHSun et al.ll201ll ). Steeper radio spectra 
may be explained by a cumulative effect of non-linear 
and oblique-shock steepening in the case of young SNRs 
( Bell. Schure fc Revilld 120111 ). An important predic- 
tion of the so called non- linear DSA (NDSA), where 
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the reaction effects of cosmic-ray particle pressure is 
taken into account, is that the spectrum of young SNRs 
should flatten at higher energies so that a "concave uj 
spectrum is formed (Reynolds 2008; Helder et al. 201^ 

There is al so a considerabl e number of Galactic SNRs 
with a < 0.5 ( Revnoldsll2011 ). Contamination with flat 
spectrum thermal emission can not be completely ruled 
out in all of the cases so it may be partially responsible 
for these lower (flat) spectral index values. 

The most of the (shell, composite and mixed- 
morphology) SNRs with a < 0.5 are in fact (evolution- 
ary) old objects expanding in a high density environ- 
ment. Because the shock velocities are much lower than 
in young SNRs, the oblique-shock effects considered by 
Bell et al. (2011) are unlikely to contribute to the flat- 
tening of the radio spectral indices of these SNRs. Also, 
in the case of older SNRs, it is not yet clear whether 
the mechanism for producing relativistic electrons is lo- 
cal acceleration at a shock front or compression of the 
existing population of Galactic background relativistic 
electrons (Leahy 2006 and references therein). 

It must be noted that the majority of members of 
the so called mixed-morphology class of SNRs have 
flat spectral indices (see Table 4 in Vink 2012). The 
mixed-morphology (thermal composite) SNRs are char- 
acterized with bright (thermal) interiors in X-rays, 



and bright (n on-th ermal) shel l-like radio morphology 
( Rho fc Petrelll998l : IVinMl2012l ). They are mainly evo- 



lutionary old SNRs expanding in a high density en- 
vironment (many of them interacting with mol ecular 
cloud s) and many of them are detected in 7-rays ( Vink 
2012 ). In Greens's catalog of Galactic SNRs ([Green 



2009f l these remnants are classified either as shell or 



composite type (they are not separated from the rest 
of SNRs). 

It should be pointed out that, from the list of the 
SNRs with flat spectrum, plerions (pulsar wind or syn- 
chrotron nebulae) are excluded since they are not of 
interest for this paper. In fact, as the radio and X-ray 
emission from these objects are powered by the pulsar 
wind, not by the supernova explosion, they should not 
be referred as SNRs (at least in a classical context). For 
pulsar wind nebulae, radio spectral indices are usually 
less than around 0.3 ( Reynolds et al.ll2012l) . 

Finally, for a completeness of this summary, it should 
be noted that among the Galactic and some Large 
Magellanic Cloud (LMC) SNRs, existence of significant 
spectral break and steepening at higher radio frequen- 
cies were also observed (Crawford et al. 2008; Xiao et al. 
2008; Bozzetto et al. 2010; de Horta et al. 2012). The 
best example is the radio continuum spectrum of Galac- 
tic SNR G180.0-1.7 (S147) for which the spectral break 
is around 1.5 GHz, with low- frequency spectral index 



around 0.3 and high-frequency spectral index around 
1.20 (Xiao et al. 2008). The compression of the local 
magnetic field and shift of the turn-over in the Galactic 
radio spectrum to higher frequencies is one of the possi- 
ble explanation of this kind of spectrum. On the other 
hand, synchrotron losses during the early phase of the 
SNR would cause a bend at a rather high frequency, 
but subsequent expansion of the SNR would shift the 
frequency toward about 1.5 GHz (Xiao et al. 2008). 
Some of the models that try to explain " concave down" 
spectra of some Galactic and extragalactic objects (in- 
cluding synchrotron losses as well as finite emission re- 
gion) are presented in the following papers: Bregman 
(1985), Biermann & Strittmatter (1987), Heavens & 
Meiseheimer (1987). 

It is very difficult to derive accurate integrated flux 
densities. The small number of data points (with asso- 
ciated errors) and the dispersion of flux densities at the 
same frequencies pose a practical problem in determin- 
ing the (integrated) radio spectral index and, generally, 
the shape of radio continuum. In that sense, it should 
be noted that present knowledge of radio spectral in- 
dices as well as of the overall shape of the Galactic 
SNR radio continuum spectra is far from being precise 
(at least for the majority of Galactic SNRs). New ob- 
servations usually lead to non-negligible changes in the 
measured spectral indices. For example, altho ugh for 
some mixed-morphology SNRs Greens's catalog ([Green 
2009f ) gives ra ther steep spectr al in dex, new q b serva, - 



tions given in iGao et al.l (J2011I) and ISun et al.l (|2011[ ) 
give slightly flatter, but still steep {a > 0.5) spec- 
tra (e.g. for G53.6-2.2, a = 0.75 -^ 0.50; G93.7-0.2, 
a = 0.65 -^ 0.52; G116.9+0.2, a = 0.61 -^ 0.57; 
G160.9+2.6, a = 0.64 -^ 0.59). In this paper are con- 
sidered only th ose SNRs for which the flux densities are 
on the scale of iBaars et al.l ( 19771 ) and have known er- 
rors < 20%. SNRs with only four or less data points 
are not discussed in this paper. 

Table 1 shows the list of several Galactic SNRs, 
with well defined radio spectra in the sense of previ- 
ously mentioned criteria, and their basic properties: 
type, integrated radio spectral index a, presence of 
low-frequency turn-over in overall spectrum (t.o.), dis- 
tance d, mean diameter D, approximate age t, associ- 
ation/interaction with molecular cloud (MC), and de- 
tection in 7-rays (7). 

The question of nature of the flat spectral indices in 
some Galactic SNRs is, in fact, still open. In Sections 
2 - 4 a brief revision of proposed explanations (with 
the strongest arguments in their favor) is presented. 
Section 4 also deals with high energy properties of the 
SNRs with flat radio spectra focusing on their detec- 
tion in 7-rays with the Large Area Telescope onboard 



Table 1 Properties of several Galactic SNRs with a < 0.5 that expand in high density environment: type, spectral index a, 
low-frequency turn-over in overall spectrum (t.o.), distance d, mean diameter D, approximate age t, association/interaction 
with molecular cloud (MC), and detection in 7 rays (7). 



SNR 


Type 


a 


t.o. 


d [kpc] 


D[pc] 


t [kyr] 


MC 


7 


Ref. 


G6.4-0.1 (W28) 


M-M 


0.35* 


7 


1.9 


26.5 


33-150 


Y 


Y 


1-5 


G18.8-H0.3 (Kes 67) 


S 


0.40 


N 


14 


57 


16-100 


Y 


? 


6-8 


G31.9-f0.0 (3C391) 


M-M 


O.49I" 


Y 


8 


-14 


4 


Y 


Y 


9-12 


G34.7-0.4 (W44) 


M-M 


0.37 


N 


2.9 


-26 


20 


Y 


Y 


13-16 


G39.2-0.3 (3C396) 


C 


0.34 


Y 


6.2 


-13 


3 


Y 


7 


13, 17-19 


G43.3-0.2 (W49B) 


M-M 


0.46 


Y 


8-11 


9-13 


2.3 


Y 


Y 


13, 20-22 


G89.0-H4.7 (HB21) 


M-M 


0.36 


N 


0.8 


-25 


19 


Y 


Y 


23-25 


G94.0-hl.O (3C434.1) 


S 


0.45 


N 


4.5 


-39 


25 


Y? 


7 


13, 26-27 


G189. 1-1-3.0 (IC443) 


M-M 


0.38 


Y 


1.5 


20 


4 


Y 


Y 


23, 28-31 



(1) Abdo at al. (2010a), (2) Dubner et al. (2000), (3) Gabici et al. (2010), (4) Giuliani et al. (2010), (5) Sawada & 
Koyama (2012), (6) Tian et al. (2007), (7) Paron et al. (2012), (8) Dubner et al. (1999), (9) Brogan et al. (2005), (10) 
Chen et al. (2004), (11) Castro & Slane (2010), (12) Green (2009), (13) Sun et al. (2011), (14) Abdo et al. (2010b), (15) 
Castelletti et al. (2007), (16) Uchida et al. (2012), (17) Su et al. (2011), (18) Patnaik et al. (1990), (19) Anderson & 
Rudnick (1993), (20) Abdo at al. (2010c), (21) Moffett & Reynolds (1994), (22) Zhou et al. (2011), (23) Gao et al. (2011), 
(24) Reichardt et al. (2012), (25) Pannuti et al. (2010), (26) Foster (2005), (27) Jeong et al. (2012), (28) Troja et 
al. (2006), (29) Troja et al. (2008), (30) Castelletti et al. (2011), (31) Abdo et al. (2010d) 



* ICreenl (|2009l ') noted that spectral index of this SNR is varying. 

t Or 0.54 at frequencies higher and 0.02 at frequencies lower than 1 GHz (Sun et al. 2011). 



Fermi 7-ray Space Telescope (Fermi LAT). In Section 
5, model which includes significant intrinsic thermal 
bremsstrahlung radiation is discussed. Finally, Section 
6 concludes the paper. 



2 Models involving second-order Fermi 
mechanism 

Although Fermi I (DSA) is generally believed to be a 
primary mechanism for particle acceleration in SNRs, 
Fermi II (stochastic acceleration, see Liu et al. 2008 for 
details) process in the shock turbulent vicinity could 
al so be significant in so r ne cas es . 

ISchlickeiser fc Fiirstl (|l989l) used results of Droge et 
al. (1987) who solved the steady-state transport equa- 
tion for the volume-integrated phase space distribution 
function of relativistic electrons in the vicinity of a 
plane-parallel shock wave allowing simultaneously for 
energy gain by multiple shock crossing and by second- 
order Fermi momentum diffusion upstream and down- 
stream of the shock. They concluded that the observed 
dispersion in spectral index values below a = 0.5 is 
attributed to a distribution of low upstream plasma /? 
values for different SNRs. 



Plasma /3 is defined as follows: 
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where P is the gas pressure, Pmag the magnetic pres- 
sure, B the magnetic induction, n the number density, 
k the Boltzmann constant, T^ and T{ the electron and 
ion temperatures (if different, or just temperature T of 
the medium), T4ound the adiabatic sound speed, Va the 
Alfven speed and 7 the adiabatic index. For example, 
temperatures of around 10^ K (— 1 eV) and number 
densities of around 0.1 cm^"^ lead to the plasma j3 of 
0.16 for the average value of Galactic magnetic field 
(magnetic induction) of around 5 /iG. Higher densities 
(as well as higher temperatures) would rise the plasma 
j3. Also, higher magnetic fields will result in lower val- 
ues of plasma (3 for the same set of temperatures and 
de nsities. 

Schlickeiser &: Fiirstl (|l989l ) also pointed out that 
values of /? ~ 0.05 are sufficient to yield a ~ 0.2 for 
shock compression ratios higher than (or equal) 2.5 (as- 
suming, also, constant spatial diffusion coefficient; see 
equations 6 and 7 in Schlickeiser & Fiirst 1989). 

On the other hand, in their work, they did not ac- 
count for the possibility of compression ratios higher 



than 4 that could be expected in some (evolutionary) 
older SNRs (spectral index becomes 0.5 for compres- 
sion ratio of 4 independently of plasma /?). Droge et 
al. (1987) emphasized that the flat spectra could be ob- 
tained for compression ratios that tend to the value of 4 
if the spatial diffusion coefficient increases with momen- 
tum and /3 = 1. However, in the case of /3 = 1 plasma 
and constant spatial diffusion coefficient, flat spectral 
indices are obtained for compression ratios higher than 
4 (which is not discussed in Droge et al. 1987). Adi- 
abatic index is assumed to be 5/3 in all the relevant 
eq uations. 

Ostrowskil ()l999r ) discussed the particle acceleration 



process at a (parallel) shock wave, in the presence of 
the second-order Fermi acceleration in the turbulent 
medium near the shock, as an alternative explanation 
for the observed flat radio spectra of some Galactic 
SNRs in molecular clouds. He found that even very 
weak shocks may produce very flat cosmic-ray particle 
spectra in the presence of momentum diffusion (using 
power-law form for the spatial diffusion coefficient - 
K oc p°'^^, assuming shock compression ratio of 4 and 
applying an effective jump condition for the Alfven ve- 
locity - Va,i = Va,2, where indices 1 and 2 refer to the 
upstream and downstream of the shock, respectively). 
Ostrowski (1999) pointed out that as a consequence, the 
dependence of the particle spectral index of the shock 
compression ratio can be weaker than that predicted 
for the case of pure first-order acceleration. His analy- 
sis was based on the model of Ostrowski & Schlickeiser 
(19 93) that improved the work of Droge et al. (1987) 
and lSchlickeiser fc Fiirsti ( 19891 ) deriving a simplified ki- 
netic equation in the momentum space valid for the case 
of small momentum diffusion and also allowing the de- 
sc ription of s t ationa ry spectra at the shock. 



Ostrowskil (|1999l ) emphasized that shock waves with 



Alfven speed non-negligible in comparison to the shock 
velocity are responsible for generating the flat parti- 
cle distribution. He noted that for the ratio of Alfven 
speed and shock velocity (inverse of the so called Alfven 
Mach number, see equation 4) around 0.1 the second- 
order Fermi process can substantially modify the energy 
spectrum of the particles accelerated at the shock and 
that the model holds for values of that ratio less than 
0.2. 

Previously mentione d Alfven M ach number of a 
shock is determined by (|Bvkovll2004l) : 



Ma 



Va 



460 148 ni/2 / B_ 



(4) 



where Vs is the shock velocity, Va is the Alfven velocity 
and 1^8 is the shock velocity in 10^ cm s~^, na is the 
ionized ambient gas number density in cm~'^ and B^q 



is the local magnetic induction, just before the shock, 
measured in fiG. Alfven Mach number is generally an 
important parameter as the ratio of acceleration rates 
of the Fermi I (DSA) a nd Fermi II (sec ond-order) pro- 



cesses scales as Ma ^ (iRevnoldsl 1201 It ) . As the bulk 



post-shock flow is expected to be super-Alfvenic it is 
logical to conclude that the DSA is more rapid process. 

The corresponding analysis was successfu Uy applied 
to the SNRs W44 and IC443. Ahhough lOstrowski 
(|l999l) could not reproduce the steep spectral index 
(0.55) in the case ofSNR 3C391, using the spectral 
index value of 0.49 (JGreenl 120091 ) . this model could be 
justified (for shock velocity of around 300 km/s and 
inter-clump density of around 10 cm~^, as in Ostrowski 
1999, this model could account for observable spectral 
index of 0.49 if magnetic induction is ~ 17 /iG; see 
equations 2.6-2.9 from Ostrowski 1999). It is worth 
mentioning that compression ratios higher than 4 (not 
discussed in Ostrowski 1999) would give, in this model, 
lower magnetic induction estimates. 

It is worth noting that SNRs expanding in molecu- 
lar clouds mostly evolve in the clumpy medium. The 
dominant part of the cloud mass is contained in the 
compact dense clumps filling only around 10% or less of 
the cloud volume and the rest of the cloud is filled with 
inter-clump gas of average number density of 10 cm~'^ 
(Chevalier 1999; Ostrowski 1999). However, one must 
bear in mind that SNR precursors could contribute to 
modification of physical conditions in the immediate 
pre-shock ISM (Arthur 2012). 

The fact that many SNRs with flat spectral indices 
are also connected with molecul ar cloud environm ent 
could be a sign that models like lOstrowskil ( 19991 ) are 
more appropriate in these cases. As the radiative shock 
has been revealed by the Spitzer IRS slit observation 
of SNR 3C396, that expands through the clumpy en- 
vironment, model of Ostrowski (1999) could possibly 
account for the observed flat radio spectral index (He- 
witt et al. 2009; Su et al. 2011). For inter-clump density 
of around 1 cm^'^ and for shock velocities of 120 km/s 
and 350 km/s (Hewitt et al. 2009), magnetic induc- 
tion of 10 — 29 /iG is needed, respectively, to account 
for observable radio spectral index of 0.34 (see equa- 
tions 2.6-2.9 from Ostrowski 1999). For a blast wave 
velocity of 870 km/s (Su et al. 2011), needed value 
of magnetic induction is higher - around 72.5 /iG. In 
the case of SNR W28, for shock velocity of around 80 
km/s and inter-clump density of around 5 cm~^ (Rho & 
Borkowski 2002), model of Ostrowski (1999) could ac- 
count for observable spectral index of 0.35 if magnetic 
induction is ~ 14 /iG. For higher inter-clump density of 
10 cm~'^, needed magnetic induction would be around 
20 /iG. Finally, HB21 is known to be interacting with 



clumpy molecular clouds (Pamiuti et al. 2010 and refer- 
ences therein). If we assume, approximately, 130 km/s 
for the velocity of radiative shell and the mean ambient 
Hi density of 3.7 cm~'^ (Koo et al. 2001 and references 
therein), model of Ostrowski (1999) could account for 
observable radio spectral index of 0.36 if magnetic in- 
duction is '^ 19 /iG. In all of these cases compression 
ratio of 4 was assumed. 

Finally, it must be emphasized that there are sev- 
eral SNRs, with flat spectral indices, that are gener- 
ally associated with relatively lower densities such as 
G82.2+5.3 (W63), with a = 0.44 (Higgs et al. 1991; 
Mavromatakis 2004; Gao et al. 2011), and G166.0-I-4.3 
(VRO 42.05.01), with a = 0.33 (Guo & Burrows 1997; 
Leahy & Tian 2005; Bocchino et al. 2009; Gao et al. 
2011). Both SNRs are classified as mixed-morphology 
ones. W63 is probably expanding in the wind blown 
bubble and VRO 42.05.01 has probably broken out 
of the cloud within which it formed, expanded across 
an interstellar tunnel or cavity, and is now interact- 
ing with the material that forms the opposite tun- 
nel wah (Pineauh et al. 1987; Burrows & Guo 1994; 
Guo & Burrows 1997; Mavr o matak is 2004). If the 
model of ISchlickeiser fc Fiirsd (|l989() is acceptable for 
these SNRs, then the observed spectral indices lead 
to the plasma /3 around 0.17 and 0.05 for W63 and 
G166.0-I-4.3, respectively, in the case of compression 
ratio of 3.9 and constant spatial diffusion coefficient. 
Slightly higher values of plasma j3 are obtained for com- 
pression ratio of 2.5. If power-law index of the spatial 
diffusion coefficient dependence of momentum equals 
0.67, then observed spectral indices could not be ob- 
tained by the model of Droge et al. (1987). Finally, 
if plasma P — \ than compression ratios of around 
6.95 and 12.3 correspond to the observed spectral in- 
dices of W63 and G166.0-I-4.3, respectively (for con- 
stant spatial diffusion coefficient). These estimates will, 
on the other hand, lead to (unrealistically) low isother- 
mal Mach number, fully radiative SNR shocks. 



3 Models involving high compression ratios 

The flatter spectra may also be due to compression ra- 
tios greater than four at fully radiative (isothermal) 
shocks (in the post Sedov- Taylor phases of SNR evo- 
lution) . In the case of parallel isothermal shock, com- 
pression ratio X is equal to a square of isothermal Mach 
number (Mt) so that test particle DSA gives: 



a 



2(x-1)' 



X 



M^. 



(5) 



only for very low isothermal Mach number fully radia- 
tive shocks (isothermal Mach numbers around 2-3) in 
the test-particle DSA. 

Of course, one must bear in mind that the onset 
of the so called momentum-conserving snowplow phase 
occur only near the end of a SNR lifetime, or not oc- 
cur at all and even a pres sure-driven snowplow mod el is 
just the asymptotic case ( Bandiera fc Petrukll2004 ). Of 
course, local regions of the SNR blast wave may leave 
the Sedov- Taylor phase earlier at those places where the 
ambient density is much higher than average, though 
the bulk e volution still cor responds to energy conserv- 
ing phase (|Revnold£ll2008l ). This means that different 



parts of an SNR may be in different phases so that dif- 
ferent shock velocities and compression ratios within 
the SNR are possible, which could have repercussions 
on the integrated radio spectrum. 




In Figure 1 it is shown that this explanation, for observ- 
able flat spectral index range, a G (0.25 — 0.5), holds 



Fig. 1 The change of integrated radio spectral index with 
isothermal Mach numbers of fully radiative (isothermal) 
SNR shock wave. 

It must be noted that the values of adiabatic index 
slightly less than 5/3 (due to high temperatures be- 
hind the shock front and/or because the molecules in 
the cloud, with which SNR interacts, will cool more 
efficiently than the ambient medium or the ejecta as 
discussed in Davis 2007) would reproduce some of the 
observable flat spectra in the case of parallel adiabatic 
shock wave and test-particle DSA (see equation 2). On 
the other hand, such values are generally not expected 
in the post-shock medium. Finally, Anand (2012) de- 
rived and discussed the jump relations across a shock in 
non-ideal gas flow. It is interesting to note that even for 
very small non-idealness of gas, using equation (5) from 
Anand (2012), for the compression ratio, flat spectral 
indices are obtained for the large range of shock veloc- 
ities in the framework of DSA. Of course, applicability 



of the non-ideal gas equation of state in the case of the 
particular ISM plasmas is questionable. 



4 High energy properties of the SNRs with flat 
integrated radio spectral indices 

Since there is a significant connection between the ma- 
jority of Galactic SNRs with flat integrated radio spec- 
trum and their detection in 7-rays, the high energy 
properties of these SNRs are now discussed (see Table 
1). Apart from historical (e.g. Tycho and Cassiopeia A 
- Abdo et al. 2010e; Giordano et al. 2012) and young 
TeV-bright SNRs (e.g. RX J1713. 7-3946 and Vela Jr. - 
Abdo et al. 2011; Tanaka et al. 2011), among the SNRs 
recently detected with Fermi LAT, there is a huge class 
of SNRs that interact with molecular clouds as well as 
of the evolved SNRs without molecular cloud interac- 
tion (e.g. Cygnus Loop, S147 - Katagiri et al. 2011; 
Katsuta et al. 2012). 

7-ray emission from SNRs can be produced by few 
different, non-exclusive, mechanisms. Electrons and 
positrons in the SNR shell can interact by inverse 
Gompton scattering with ambient photon fields (such 
as cosmic microwave background and infrared radia- 
tion) to produce 7-rays (Helder et al. 2012). Also, 
electron-atom or electron-molecule interactions in a 
dense medium result in bremsstrahlung radiation (Re- 
ichardt et al. 2012). Mentioned mechanisms involve 
electrons accelerated up to the sufficient energies to pro- 
duce the observed 7-rays (leptonic scenario). 7-rays can 
also be produced by neutral pion decay, where neutral 
pions result from the collision of accelerated protons 
(or heavier nuclei) with nucleons of the ambient gas 
(hadronic scenario). 

In the case of SNRs that interact with molecular 
clouds, hadronic scenario is probable and one of the 
spectral signatures is steepening in the GeV band. Gen- 
erally, absence of non-thermal X-ray emission favors 
hadronic origin of the observed 7-ray emission. The 
rapid steepening of the spectrum above few GeV is the 
kind of signature expected from the re-acceleration of 
pre-existing cosmic-rays, where high energy cosmic-rays 
escape from the SNR confinement region (Reichardt et 
al. 2012 and references therein). Uchiyama et al. (2012) 
emphasized that it has not been clear yet whether the 
7-ray emission is produced by escaping or by trapped 
cosmic-rays, since shocked molecular clouds inside an 
SNR could also be the sites of efficient 7-ray produc- 
tion. Although the majority of SNRs that interact with 
molecular clouds and are detected in 7-rays are those 
with flat radio spectra (e.g. W44, IC443, W28, W49B, 
3G391, W51C, GTB37A, HB21 - Abdo et al. 2009; 



Abdo et al. 2010abcd; Gastro & Slane 2010; Reichardt 
et al. 2012), it must be mentioned that several SNRs 
that interact with molecular clouds and have spectral 
indices of 0.5 are also detected in 7-rays with Fermi 
LAT (e.g. G8. 7-0.1, G109. 1-1.0, G349.7+0.2 - Gastro 
& Slane 2010; Ajello et al. 2012; Gastro et al. 2012). 

Reichardt et al. (2012) emphasized that SNR HB21 
belongs to the group of low-luminosity, GeV-emitting 
SNRs (such as Gygnus Loop or S147), which are clearly 
less luminous than the first GeV-emitting SNRs that 
were discovered (such as W51G, IG443, W49B). The 
spectral break is found at lower energy in HB21 than 
in the case of the luminous SNRs. Reichardt et 
al. (2012) suggested that the 7-ray emission from IIB21 
can be understood as a combination of emission from 
shocked/illuminated molecular clouds, one of them co- 
incident with the SNR shell itself. 

SNR G260.4-3.4 (Puppis A), with radio spectral in- 
dex of 0.5 (Green 2009), was also observed at 7-rays 
(Hewitt et al. 2012) and it represents an interesting 
transitional case between young SNRs still evolving into 
a circumstellar medium (e.g. Gas A), and older SNRs 
which are interacting with large, dense molecular clouds 
(e.g. IG 443). 

In the case of Kes 67, not yet detected by Fermi 
LAT, the bright pulsar wind nebula HESS J1825-137 
is extended near the position of the SNR (Grondin et 
al. 2011) thus making the analysis highly complicated. 
Similarly, near the position of SNR 3G396, a radio-faint 
7-ray pulsar (PSR J1907-h0602) powering a bright TeV 
pulsar wind nebula (MGRO J1908-t-06) is located, sub- 
stantially complicating the analysis (Abdo et al. 20101). 

Possible interaction with molecular cloud in the case 
of SNR 3G434.1 is indicated by Jeong et al. (2012), 
which makes this remnant an interesting target for anal- 
ysis of possible 7-radiation in its direction. On the other 
hand, there are no point sources included in the LAT 
2-year point source catalog (Nolan et al. 2012) closer 
than 2 degrees to the position of SNR 3G434.1. Also, 
there are no extended sources closer than 10 degrees to 
the position of this remnant. A detailed study of the 
7-emission in the direction of SNR 3G434.1 is planned 
for a work in the near future. 

4.1 Model of Uchyama et al. 



Uchiyama et al.l (J2010|) analyzed 7-ray emission from 
SNRs interacting with molecular clouds. They showed 
that the simple re-acceleration of pre-existing cosmic- 
rays by the process of diffusive shock acceleration at 
a cloud shock is generally sufficient to power the ob- 
served 7-ray emission through the decays of neutral 
pions. Neutral pions are produced in hadronic in- 
teractions between high-energy protons (nuclei) and 



gas in the compressed-cloud layer. lUchivama et al, 
( 20101 ) proposed that the radio emission may be ad- 
ditionally enhanced by the presence of secondary elec- 
trons/positrons, i.e. the products left over from the de- 
cay of charged pions, created due to cosmic-ray nuclei 
colliding with the background plasma. They concluded 
that presence of secondary electrons/positrons may also 
explain the flat spectral radio indices of some mixed- 
m orphology SNRs. 



Uchivama et al. 



( 20101) applied their analysis to 
SNRs W51C, W44 and IC443. Of course, this kind 
of explanation holds only for those SNRs that interact 
with molecular clouds assuming further that hadron 
scenario for 7-radiation is significant. Another restric 



tionof this model is the fact that lUchivama et al. 



( 20101 ) used relations that are applicable in the case 
of the Sedov- Taylor phase which are not suitable for 
mixed-morphology SNRs. On the other hand, the sim- 
ple re-acceleration of pre-existing cosmic-rays and sub- 
sequent compression alone would not fully explain the 
7-rays associated with cloud-interacting SNRs. It must 
be noted that GeV and TeV 7-rays outside the southern 
boundary of SNR W28 may be explained by molecular 
clouds illuminated by runaway cosmic -rays (Uchiyama 
et al. 2010 and references therein). lUchivama et al. 



( 20101 ) also assumed pre-existing cosmic-rays in the 
cloud to have the same spectra as the Galactic cosmic- 
rays in the vicinity of the Solar system although the 
ambient cosmic-rays in the pre-shock cloud may deviate 
from the Galactic pool due to the runaway cosmic-rays 
that have escaped fr om SNR shoc k s at e arlier epochs. 
Another model bv lBvkov et al.l ( 20001 ) was also pro- 
posed assuming the possibility of electron acceleration 
by magnetohydrodynamical turbulence behind high 
density shocks (Vink 2012). In their model, the calcu- 
lated radio spectrum fits that observed from the SNR 
IC443 shell because of the large shock compression ra- 
tio and to a lesser extent because of the effect of the 
second-order Fermi acceleration. It is worth mention- 
ing that this model involves a lepton scenario of 7-ray 
emission from SNRs interacting with molecular clouds. 



5 Model involving significant intrinsic thermal 
bremsstrahlung radiation 

Recently, it was shown that the observations over a 
very broad range of radio frequencies reveal a curva- 
ture in the spectra of some evolutionary older Galactic 
SNRs expanding in high density environment (Tian & 
Leahy 2005; Urosevic & Pannuti 2005; Leahy & Tian 
2006; Urosevic, Pannuti & Leahy 2007; Onic & Urosevic 
2008). The NDSA effects could not be responsible 



for a characteristic shape of the integrated radio spec- 
tra in these cases. The presence of intrinsic thermal 
bremsstrahlung radiation was proposed as an explana- 
tion of the "concave up" radio spectrum (Urosevic 2000; 
Urosevic et al. 2003ab; Onic et al. 2012 and references 
therein). Presence of the cooled thermal X-ray elec- 
trons in post Sedov- Taylor SNRs, higher ambient den- 
sities, interaction with or expansion in molecular cloud 
are main arguments used by the model that includes 
significant intrinsic thermal bremsstrahlung emission 
(thermal model). Detections of the low- frequency turn- 
overs, in the integrated radio spectrum, related with 
thermal absorption linked to SNRs (Brogan et al. 2005; 
Castelletti et al. 2011) are important in justification 
of this model as one can then predict the significance 
of thermal emission at higher frequencies. Generally, 
linear polarization measurements (that give the upper 
limits for the thermal component of the radio emission), 
possible detections in Ha as well as the observations of 
radio recombination lines associated with SNRs should 
accompany this kind of analysis. 

It must be emphasized that most of the potential 
targets for testing the thermal model hypothesis fall 
in the mixed-morphology class of SNRs. Because of 
the fact that these SNRs appear to be in the radiative 
phases of their evolution, with shock velocities less than 
200 km/s (Vink 2012), the ensemble of cooled ther- 
mal X- ray electrons coul d exist, supporting the thermal 
model ( Onic et al.ll2012 i). In addition, recent discovery 
of strong radiative recombination continua (RRC), seen 
in the X-rays of several SNRs of this type, represents 
definite evidence that their plasma is recombining (Ya- 
maguchi et al. 2012). RRCs in their X-ray spectra are 
indicative of higher densities in interiors of these rem- 
nants (Vink 2012). Thus, mixed-morphology SNRs are 
characterized by a more/less uniform high interior den- 
sity with medium hot temperatures, rather than very 
high temperatures and very low interior densities, as 
in the case of typical shell remnants. Also, among the 
SNRs with detected RRC, the most of them (if not 
all) have fiat radio spectral index (Sawada & Koyama 
2012; Uchida et al. 2012; Vink 2012). Finally, as it 
was emphasized earlier, the morphology of these SNRs 
is diffic ult to expla in with standard SNR evolutionary 
phases ( Vinkll2012t ) and several models were proposed 
in the current literature (White & Long 1991; Cox et 
al. 1999). An X-ray emitting overionized plasma is the 
result of an early heating followed by a rapid cooling, 
but the mechanism responsible for such processes in 
mixed-morphology SNRs is still not completely under- 
st ood (Miceli 2011) 



Onic et al.l (|2012|) have recently analyzed the inte- 



grated radio spectra of 3 SNRs (IC443, 3C391, 3C396) 



for which the thermal model could be a natural expla- 
nation. All of the 3 SNRs are characterized by flat spec- 
tral indices (see Table 1) based on a simple power-law 
(synchrotron) fit and all of them interact with molec- 
ular clouds. For IC443 and 3C39f the low-frequency 
turn-over is due to thermal absorption linked to the 
SNRs and for SNR 3C396 the nature of the observed 
low-frequency turn-over is not yet clear. It must be 
emphasized that most of the SNRs that are the best 
targets for testing thermal model are also characterized 
by the flat spectral indices. The observed flat spectral 
indices could be just the apparent manifestation of the 
presence of thermal component. It must be noted that 
there are no predictions of curved - "concave up" in- 
tegrated radio spectra in the cases of the explanations 
presented in Sections 2-4 for the flat radio spectra seen 
in some SNRs. 

The shock wave dissociates molecules, or raises the 
gas temperature so that previously inaccessible degrees 
of freedom become accessible. The compression ratio 
for the parallel shock wave, including the possibility of 
jump in adiabatic index, is as follows: 
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where 71 and 72 are adiabatic indices in upstream and 
downstream region, respectively. If 71 = ^ (diatomic 
molecules like II2) and 72 = I (fully ionized gas) it 
is easily seen that (upstream) adiabatic Mach num- 
ber Af = 20 leads to a sa 0.51 so that for weaker 
shocks the steeper spectra will be preferable. As the 
ambient gas is usually pre-ionized by radiative precur- 
sor, the similar conclusion (for slightly smaller adia- 
batic Mach numbers) could be drawn from the sim- 
ple relation for constant adiabatic index of 5/3 (de- 
scribed in Section 1 using equation 2, see Figure 2). 
This means that we expect that evolutionary old SNRs, 
with velocities smaller than around 200 km s""'^, have 
steeper synchrotron spectral indices in the framework 
of test-particle DSA. This can partially account (except 
high parameter uncertainties) for steeper synchrotron 
spectral indices det ermined from the thermal plus non- 
thermal model fits ( Onic et al.ll2012 ). 

Radio spectral index variations, both with frequency 
and location within several Galactic SNRs were de- 
tected (Fiirst & Reich 1988; Anderson & Rudnic 1993; 
Zhang et al. 1997; Leahy & Roger 1998; Leahy et al. 
1998; Tian & Leahy 2005; Leahy 2006; Ladouceur & 




Fig. 2 Change of the radio spectral index with Mach num- 
ber for slow shocks using equations (2) - full line and (6) - 
dotted line. 



Pineault 2008). In order to have spectral changes with 
position, physical conditions that alter the radio spec- 
trum must change with position. Of course, spatial ra- 
dio spectral index variations are naturally expected in 
the case of composite SNRs due to the presence of pul- 
sar wind nebula. On the other hand, variations in a are 
detected in both shell and mixed-morphology SNRs. If 
the spectral index is changing with frequency, generally, 
there is a possibility that two different radiation mecha- 
nisms are producing the radio continuum emission. In- 
tegrated radio spectrum with " concave up" curvature in 
the case of older SNRs, where efficient NDSA process is 
not expected, can be produced by adding two different 
emission spectra along the same line of sight or within 
the same emitting volume (Leahy & Roger 1998). These 
kind of spectra could be due to different electron pop- 
ulations or because of the significant intri nsic thermal 
bremsstrahlung radiation ( Onic et al.l2012 ). Linear po- 
larization measurements are of the great importance as 
spectral flattening associated with high density regions 
of an SNR which do not show high linear polarization 
could be naturally explained by the thermal model. 

It is worth mentioning that there are some other 
processes that could shape the integrated radio con- 
tinuum at higher frequencies (between 10-100 GHz) in 
several Galactic SNRs, such as the spinning dust emis- 
sion (Draine & Lazarian 1998; Scaife at al. 2007). The 
spinning dust model will be responsible for a creation 
of a "hump" in the high-frequency radio spectrum (see 
Figure 5 from Scaife at al. 2007). Of course, as it was 
emphasized it the Section 1, present knowledge on in- 
tegrated radio spectrum is generally unsatisfactory for 
firm quantitative analysis of the proposed models so 
as to the contribution of dust emission. More data at 



radio frequencies higher than 1 GHz, e.g. those from 
ATCA (The Austraha Telescope Compact Array), are 
expected to shed a hght on the existence of the so called 
"radio thermally active" SNRs as well as on the issue 
of the flat spectral indices. It must be mentioned that 
there is a general problem in observing at the higher 
radio continuum frequencies (>30 GHz) due to atmo- 
spheric absorption issues. Ground-based radio astron- 
omy is limited to high altitude sites such as in the case 
of ALMA ( Atacama Large Millimeter Array) . 

5.1 The results regarding thermal model for Galactic 
SNRs 

Except the supernova remnants G31. 9-1-0.0 (3C391), 
G39.2-0.3 (3 C396 ), G189.1-F3.0 (IC443), analyzed in 
Onic et al.l (|2012l) . that represent the best targets 
for justification of the thermal model, in this pa- 
per, the integrated radio spectra of SNRs G6. 4-0.1 
(W28), G18.8-I-0.3 (Kes 67), G34.7-0.4 (W44), G43.3- 
0.2 (W49B), G89.0-H4.7 (HB21), G94.0+1.0 (3C434.1) 
were also discussed. All of these SNRs are associated 
with dense environment (see Table 1). Integrated flux 
densities are collected from the literature keeping in 
mind the criteria defined in Section 1. In the case of 
presence of spectral turn-over at low frequencies only 
data which are not influenced by the possible thermal 
absorption (or maybe by sync hrotron self-absor ption) 
were used in thermal model fits ( Onic et al.ll2012l ). Sim- 
ple sum of non-thermal and thermal components, rep- 
resented by corresponding power-laws, (see equation 18 
and corresponding discussion in Onic et al. 2012) is as- 
sumed. For frequencies in GHz, the relation for the 
integrated flux density can be written as follows: 
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where: S'j^qhz ^^^ ^fcHz ^^^ ^^^ densities correspond- 
ing to thermal and non-thermal components, respec- 
tively. Weighted least squares fit (infiuenced by the 
instrumental errors) were applied. In the case of all 
of these six SNRs thermal model fit is either not sig- 
nificant (contribution of the thermal component at 1 
GHz is near zero) or the fit quality is not satisfactory 
(and associated parameter errors exceed values of the 
corespondent parameters). Pure non-thermal fit pa- 
rameters for these SNRs correspond to ones given in 
the present literature (see Table 1). 

In the case of SNR W28, low- frequency spectral turn- 
over in the integrated spectrum is possible as flux den- 
sity at 30.9 MHz has a rather low value in the sense of 
scatter from the power-law fit (see Figure 5 in Dubner 
et al. 2000). The associated error for that flux density is 



20% (Kassim 1989). On the other hand, there is a great 
spread in fiux density values, especially at 2.7 GHz, 
which prevents firm conclusions. The great spread in 
the integrated flux densities at the same frequencies is 
very good seen in the case of SNRs W49B and W44 
(Morsi & Reich 1987; Kassim 1989; Taylor et al. 1992; 
Kovalenko et al. 1994; Lacey et al. 2001; Castelletti et 
al. 2007; Sun et al. 2011). For SNR W44 there is no 
apparent low-frequency spectral turn-over in the inte- 
grated radio spectrum, only the localized one, towards 
the southeast border of the SNR, likely due to the free- 
free absorption from ionized gas in the post-shock re- 
gion at the SNR/molecular cloud interface (Castelletti 
et al. 2007). Detected turn-over at lower frequencies for 
W49B is likely due to the extrinsic free-free absorption 
by an intervening cloud of thermal electrons (Lacey et 
al. 2001). There is no apparent low-frequency spec- 
tral turn-over in integrated spectra of SNRs Kes 67, 
HB21 and 3C434.1 (Reich et al. 1983; Goss et al. 1984; 
Landecker et al. 1985; Kassim 1989; Milne et al. 1989; 
Tatematsu et al. 1990; Kovalenko et al. 1994; Dubner 
et al. 1996; Zhang et al. 2001; Reich et al. 2003; Fos- 
ter 2005; Kothes et al. 2006; Gao et al. 2011; Sun et 
al 2011). 

In the case of other SNRs that are listed in ICreenl 
( 20091) as those with flat spectral indices (not plerions), 
expanding in high density environment, there are not 
enough data (under criteria presented in Section 1) for 
a detailed discussion. Ge nerally, there is a signiflcant 
number of SNRs listed in ICreenl (|2009l ) for which the 
spectral indices are determined only by two or three 
data points. Some of the SNRs that could be interest- 
ing targets for analysis, in future, when more reliable 
data at different frequencies are obtained, are G49.2-0.7 
(W51C) and G290.1-0.8 (MSH 11-61A). 

A few words of caution regarding the justification of 
thermal model. SNR's expansion in complex environ- 
ment and possible association with H ii regions (and so 
possible thermal contamination) makes the study more 
complicated. The best example is the case of steep 
radio spectrum SNR G132.7-hl.3 (HB3) for which the 
intrinsic thermal radiation was dismissed due to the 
actual overlap between the SNR and corespondent H ii 
regions (Shi et al. 2008). Of course, the question on 
the presence of significant thermal radio emission from 
SNR HB3, in fact, still remains open as the possible 
H II region radiation contamination can not rule out the 
possibility of thermal bremsstrahlung emission from the 
SNR itself, it can just mask it. 

Final conclusion is that our present knowledge on 
integrated flux densities prevents us from possibility 
of strong justification of the thermal model and lim- 
its the_discussion_to just a few Galactic SNRs analyzed 
in lOnic et all (|2012l ). 
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6 Conclusions 

A considerable fraction of Galactic SNRs are charac- 
terized by the flat spectral indices (a < 0.5). In this 
paper, known models are summarized and discussed: 

1 . There are several explanations of the observable flat 
radio spectra of some SNRs. Most of the mod- 
els involve a significant contribution of second-order 
Fermi mechanism, but some of them discuss higher 
shock compressions, contribution of secondary elec- 
trons/positrons left over from the decay of charged 
pions, as well as the possibility of thermal contami- 
nation. 

2. In the case of expansion in high density environment, 
thermal bremsstrahlung could theoretically shape 
the spectrum of SNRs. Lack of more high quality 
data constraints the discussion only to several SNRs. 
On the other hand, this model can naturally account 
to observable curved - " concave up" radio spectra of 
some Galactic SNRs. New observations are expected 
to create a clear picture about the high frequency 
range (1-100 GHz) of the radio continuum, as well 
as about the question of the significant contribution 
of intrinsic thermal bremsstrahlung radiation. With 
instruments like the JVLA0 and ATCA, it should be 
far easier to measure the spectral indices in individ- 
ual observations. 

3. The fundamental question about the origin of fiat 
spectral indices still remains open. The more realis- 
tic picture would involve the action of more than one 
of the mentioned processes, in which some of them 
could be more or less prominent depending on the 
particular setup. 

4. Since there is a significant connection between the 
majority of Galactic SNRs with fiat integrated ra- 
dio spectra and their detection in 7-rays, as well as 
detection of RRC in their X-ray spectra, the analy- 
sis of high energy properties of these SNRs will be 
potentially very important. 
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